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(57) Abstract: A channel estimation method and apparatus for GMS/EDGE digital communications systems utilizing previously 
unexploited properties of GMS/EDGE training sequences to permit a more efficient initial channel estimation for equalizer operation 
and for joint synchronization and equalizer window sizing. In particular, any consecutive 16-symbol segment of the 26-symbol 
GMS/EDGE training sequences is both shift invariant and order invariant; and these properties enable channel estimation to be 
carried out on delayed (shifted) training sequence segments, permitting 1SI corrupted leading symbols to be avoided in computations 
and leading taps to be estimated, using the same training sequence segments, regardless of equalizer window size; and to enable all 
the 1-8 tap channels to be estimated without matrix inversion, permitting a significant reduction in computational complexity. A two- 
dimensional recursive relation is established that leads to efficient channel tap calculation for different sync points and tappositions 
which enables practical joint optimization of synchronization and window sizing for the equalizer. 
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CHANNEL ESTIMATION-BASED EQUALIZER WITH VARIABLE WINDOW LENGTH 



BACKGROUND OF THE INVENTTON 
5 1. Field of the Invention 

The present invention relates generally to the field of digital communications; and, 
more particularly, to an efficient channel estimation method and apparatus for a 
GSM/EDGE digital communications system which utilizes special properties of 
GSM/EDGE training sequences and leads also to a practical joint optimization of 
1 0 synchronization and equalizer window sizing. 
2. Description of the Prior Art 

Inter-symbol interference (ISI) is an important problem in digital communications 
systems, including those systems which operate in accordance with the Global System for 
Mobile Communications (GSM) specifications. ISI is caused by time dispersion in the 

15 transmission channel over which a signal is transmitted, and adversely affects the quality 
of the received signal. In effect, ISI causes distortion of the transmitted signal which, in 
turn, causes symbol errors in the received signal such that it becomes difficult for the 
receiver to determine what data was actually sent 

As is well-known in the art, the usual way to compensate for ISI in a GSM system 

20 is to provide a channel estimation-based equalizer in the receiver. Basically, a model or 
estimate of the propagation channel over which a received signal was transmitted is 
created; and the equalizer then uses that information to estimate the sending symbols that 
best fit the received signal. 

EDGE (Enhanced Data rates for Global Evolution) is an interface mode which has 

25 recently been developed for GSM Networks. EDGE'S principal features include new 
modulation and coding schemes which increase data capacity and speed in the air 
interface. EDGE is fully based on GSM and uses the same TDMA (Time Division 
Multiple Access) frame structure as GSM, such that it allows GSM operators to use 
existing GSM radio bands to offer wireless multimedia-based services and applications. 

30 In GSM/EDGE systems, the performance of the equalizer in combating ISI 

depends heavily on the quality of the channel estimation; and the quality of the channel 
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estimation is, in turn, highly sensitive to the accuracy of burst synchronization (the term 
"synchronization" in this document signifies "burst synchronization") and the size of the 
equalizer window. 

To effectively combat ISI, the span of the equalizer window must be large enough 
5 to cover the maximum delay spread of the channel. However, an over-sized equalizer 
window will result in performance losses in channels with shorter delay spread due to 
inaccurate channel estimation, inadequate synchronization and increased noise 
contribution. In addition, synchronization can be optimized only with the knowledge of 
the equalizer window to capture maximum energy in the window span. In general, the 

10 interdependency of synchronization and equalizer window sizing makes efficient 
optimization difficult to achieve under different channel conditions. Joint optimization was 
deemed as too expensive in practical implementation. 

Because of this interdependency, a careful compromise has to be made for a 
traditional equalizer with fixed window span so that the size of the window can provide 

15 adequate performance in long delay spread channels, such as Hilly Terrain (HT), without 
losing too much in short delay spread channels, such as Typical Urban (TU). 

When the equalizer window size is set, there are two approaches for burst 
synchronization, a correlation-based approach and a least square error (LSE) approach. 
As will be described more fully hereinafter, in current GSM receivers, synchronization is 

20 done by a correlation-based approach in which the central 16 symbols of a known 26- 
symbol training sequence is correlated with the central 16 symbols of the training 
sequence in the received signal. The synchronization is determined by either the centers of 
gravity of the correlation or the maximum correlation energy in the equalizer window. In 
the LSE approach, for each possible synchronization point, an estimation of the channel is 

25 made, and synchronization is determined by comparing the expected and received training 
sequence with least square error criterion. 

Neither of these approaches, however, is fully satisfactory. The correlation-based 
algorithm suffers from performance degradation due to inaccurate synchronization, 
especially in long dispersive channels, while a straightforward LSE-based algorithm 

30 suffers from a high degree of computational complexity (mainly due to multiple channel 
estimation). 
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SUMMARY OF THE INVENTION 

The present invention provides an initial channel estimation method and apparatus 
5 in a GSM/EDGE digital communication systems which utilizes special properties of 
GSM/EDGE training sequences to achieve a more efficient channel estimation. A two- 
dimensional recursion is further established to enable practical joint LSE optimization of 
burst synchronization and adaptive equalizer window sizing. 

More particularly, an embodiment of the present invention provides a method for 

10 efficient estimation of the transmission channel in a digital communications system which 
operates in accordance with GSM/EDGE specifications. The system has a receiver which 
receives a signal transmitted over the transmission channel and which includes a channel 
estimation-based equalizer that uses an estimate of the transmission channel to equalize 
the received signal. A received signal burst includes a 26-symbol training sequence. The 

15 method includes using any consecutive 16-symbol segment of the 26-symbol training 
sequence to estimate, without matrix operation, the transmission channel at different 
possible synchronization points for different equalizer window sizes. 

The present invention provides a more efficient channel estimation method by 
recognizing and utilizing two previously unexploited properties of GSM/EDGE training 

20 sequences which arise from the cyclic prefix structure of the training sequences. In 
particular, the present invention recognizes that the 26-symbol GSM^EDGE training 
sequences are, within certain ranges, both shift invariant and order invariant. The property 
of being shift invariant enables a channel estimation to be carried out with delayed 
(shifted) training sequence segments; which, in turn, permits: I. The ISI corrupted leading 

25 training sequence symbols to be avoided as much as possible in long dispersive channels; 
and n. The leading taps of a channel to be estimated using the same training sequence 
segments, regardless of the size of the equalizer window. The property of being order 
invariant permits channels with different time dispersion, from 1-8 symbols (as complex 
polynomials of 1-8* order) to be estimated without matrix inversion if any consecutive 16- 

30 symbol segment of the training sequence is used. This results in a significant reduction in 
computational complexity in channel estimation for varying window size which can be 
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adaptive to the actual delay spread of the channel. Combining these two properties, a two- 
dimensional recursive relation can be established for different synchronization points and 
tap positions, which makes a practical LSE joint optimization of synchronization and 
equalizer window sizing possible. 

In general, the present invention provides an efficient channel estimation for 
GSM/EDGE systems which, beside direct utilization, is fundamentally important to 
achieving a practical implementation of joint LSE synchronization and equalizer window 
sizing. 

Yet further advantages and specific details of the present invention will become 
apparent hereinafter in conjunction with the following detailed description of presently 
preferred embodiments. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 schematically illustrates the structure of a TDMA frame used in a digital 
communications system according to the GSM/EDGE standard; 

Fig. 2 schematically illustrates a "cyclic prefix" structure of GSM/EDGE training 
sequences; 

Fig. 3 is a block diagram of a simplified baseband model for a transmitter-receiver 
chain to illustrate data recovery signal processing in GSM/EDGE digital communications 
systems; 

Fig. 4 is a flow chart illustrating channel equalizer processing in a GSM/EDGE 
digital communications system; and 

Fig. 5 illustrates a two-dimensional recursive relation in channel estimation for 
different synchronization points and channel tap positions according to a presently 
preferred, embodiment of the invention. 

DETAILED DESCRIPTION OF PRESENTLY PREFERRED EMBODIMENTS 

Fig. 1 schematically illustrates the structure of a TDMA frame 20 used in a mobile 
telecommunications system according to the GSM standard. Frame 20 has a duration of 
4.615 ms and accommodates 8 information channels (time slots) 21, designated 0-7 in Fig. 
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1. Each of the 8 time slots has a duration of 0.577 ms and contains a 148-bit signal portion 
30 and an 8.25-bit guard portion (not shown) which functions to maintain separation 
between signals in adjacent time slots. The 148-bit signal portion is generally referred to 
as a "normal burst" and comprises a first 3-bit tail bit section 22, a first 57-bit coded data 
5 section 23, a first 1-bit housekeeping bit section 24, a 26-bit training sequence section 25, 
a second 1-bit housekeeping bit section 26, a second 57-bit coded data section 27 and a 
second 3-bit tail bit section 28. 

The 26-bits comprising the training sequence section 25 in a GSM digital 
communications system, is typically viewed as being divided into a centrally located 
1 0 portion 3 1 of 16-bits (sometimes referred to as the midamble portion), and side portions 32 
of 5-bits each. Alternatively, in a GSM/EDGE digital communications system, the 26-bit 
training sequence 25 can be also viewed as including a 16-bit "white" sequence 37 and a 
10-bit cyclic prefix 38 as is illustrated in Fig. 2. 



which a signal is transmitted; and in GSM/EDGE systems, the usual way to compensate 
for ISI is to provide a channel estimation based equalizer in the receiver. 

Fig. 3 is a block diagram of a simplified baseband model of a transmitter-receiver 
chain, generally designated by reference number 40, to illustrate data recovery signal 

20 processing in a GSMZEDGE digital communications system. In Fig. 3, the transmitter side 
is shown at 41 and the receiver side is illustrated at 42. A signal is transmitted from the 
transmitter side 41 to the receiver side 42 over a transmission channel 43. The 
transmission characteristics of the transmission channel frequently change due to a variety 
of factors, including the motion of the mobile terminals, the fluctuation of multi-path (time 

25 dispersive) propagation channels and variant interferences, which are generally designated 
in Fig. 3 by the noise component 44 introduced into the transmission channel. As is well- 
known to those skilled in the art, the multi-path channel and the noise component can 
adversely affect the quality of the received signal, e.g., cause intersymbol interference; and 
this necessitates that the received signal be corrected, typically by means of a channel 

30 estimation-based equalizer. 



As mentioned previously, inter-symbol interference can occur in digital 
15 communications systems as a result of time dispersion in the transmission channel over 
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Specifically, as shown in Fig. 3, at the receiver side, the received signal, which 
includes the transmitted signal with intersymbol interference and the noise component, is, 
after being filtered by filter 45, sent to an equalizer 46 which attempts to correct the 
received signal utilizing an estimate of the transmission channel as shown at 47. The 
channel estimation, in turn, is a function of burst synchronization and equalizer window 
sizing as shown at 48. The performance of the equalizer can be improved significantly it, 
in an advanced receiver, the equalizer window size is adaptive to the actual delay spread 
of the channel. Data representing an estimate of the actual transmitted symbols, as 
determined by the equalizer, is then output for further processing by the receiver. 

Fig. 4 is a flow chart which illustrates, in somewhat greater detail, the steps of the 
channel equalization processing which takes place in an advanced receiver 42 of Fig. 3. 
Specifically, after filtering the received signal at step 50, synchronization and equalizer 
window sizing according to the delay spread of the actual channel is performed as shown 
at step 52. An estimate of the transmission channel is then made at step 54, and a 
hypothetical signal provided in step 60 is applied to the estimated channel at step 56. An 
error comparison is made between the hypothetical signal applied to the estimated channel 
from step 56 and the synchronized received signal at step 58. If the result of the error 
comparison is not at a minimum value (not min output of step 58) a revised signal 
hypothesis is made in step 60, the new revised signal is applied to the hypothetical channel 
at step 56, and a new comparison is made at step 58. The comparison process is continued 
until the result of the comparison is at a minimum value (min output of step 58), at which 
time the process is completed. This processing is typically carried out with a specific 
algorithm for sequence estimation, for example, the well-known Viterbi algorithm. 

As mentioned above, and as should be apparent from Figs. 3, the performance of 
the equalizer 46 is dependent on the quality of the transmission channel estimation made 
at 47. This estimation, in turn, is highly dependent on the accuracy of the burst 
synchronization and equalizer window sizing at 48. The burst synchronization and 
equalizer window size are interdependent on one another, and this interdependency makes 
efficient optimization not trivial in different channel conditions. In contrast to the 
traditional equalizer with a fixed span, a varying window sizing adaptive to the actual 
delay spread of the channel is made possible by the efficient channel estimation proposed 
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here so that a much improved performance can be reached in both long delay spread 
channels such as Hilly Terrain, and short delay spread channels such as Typical Urban. 

Assuming an equalizer window size m, there are two approaches for burst 
synchronization; a correlation-based approach and a least square error (LSE) approach. 
5 Synchronization in current GSM receivers is done using a correlation-based 

approach by correlating the central 16 bits of the known 26-bit training sequence with the 
centrall6 bits of the training sequence of the received burst. 
20 

c(k)= Ir(i+k)t(i), k=0,...,w-l 

where w is the size of the synchronization window, r(i) is the received samples and t(i) is 
10 the known training sequence. Synchronization can be determined by either center of 
gravities of or maximum energy in the equalizer span of c(k). 

In the LSE approach, for each possible synchronization point, a least square 
estimation of the channel is made, 
H = <I>-*T T R = (T T T)- l T T R 

(2) 

15 where R is the vector of received symbols, and ^ is an mxm auto-correlation matrix of 
the training sequence symbols, 

<D = 7 r r (3) 



20 



T is a Toeplitz matrix of the training symbols 



7 = 



(4) 



Synchronization is determined by least square error criteria 

p = argmm|a"2h A )"^ ] 
where a > 1 is a penalty factor, f p (k) is an estimate of the received symbols 



(5) 
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(6) 

Neither of the above approaches, however, is fully satisfactory. The correlation- 
based algorithm suffers from performance degradation due to inaccurate synchronization. 
Performance losses are especially significant in long dispersive channels. The LSE-based 
algorithm provides higher quality channel estimation when accurate synchronization is 
obtained with equation (5); however, it has not been widely implemented due to high 
computational complexity, since, in general, equation (6) requires knowledge of all the 
channel taps at every possible synchronization position with sufficient accuracy. The 
complexity is even higher if channel estimation has to be adapted for different equalizer 
window size m using equation (2). 

In a traditional GSM receiver, only 5-tap channel estimation is usually made, 
regardless of the actual channel conditions. It is known that if the 16-symboI mid-amble 
of the training sequence is used for a 5-tap channel estimation, 
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15 7(5,1) = 

the auto-correlation matrix of the training sequences have a special form 
<l>(5,l) = r(5,l)r(5,l) = 16/ J 



(7) 



(8) 



20 



where Is is a (5 x 5) identity matrix. Here we have made the channel-taps (the order of 
channel vector polynomial), which determine the row dimension of T, and the index of the 
first training symbol in T explicit parameters of the function. Equation (8) has been used 
in S-tap channel estimation, 
1 



H = 0- , T T R=—T T R 
16 



(9) 



WO 02/103972 




12/06653 



In general, however, when Equation (2) is used for initial channel estimation for an 
equalizer window other than 5, matrix inversions is required. 

We have now discovered that the GSMZEDGE training sequences have two 
additional properties, not previously exploited, which permit channel estimation to be 
5 carried out more efficiently. More particularly, rewriting equation (8) in the following 
more general form reveals that the autocorrelation matrix of any consecutive 16-symbol 
segment of the GSM/EDGE training sequence is, within certain range, both shift invariant 
and order invariant. 
Property L Shift Invariant: 
10 ^m,d)=7*(m>d)T(m,d) = 16T OT 0<dZ(26-(l6+m)+l) (10) 

Property IL Order Invariant: 

(m y d)=T T {m,d)T(m 9 d)=\a m \<>m<.l (11) 

15 Property I is due to the cyclic prefix nature of GSM/EDGE training sequences, in which 
the 26-bit training sequence is constructed by prepending the last 10 bits to its front as 
shown in Fig. 2. This property enables a channel estimation of the form of equation (9) be 
carried out on delayed (shifted) training sequence segments. This flexible delay brings 
two important benefits: 
20 (1) ISI corrupted leading symbols can be avoided in the computation; and 

(2) the leading taps of a channel can be estimated using the same training sequence 
segments, regardless of the size of the equalizer window. 

Property II enables channel estimation to be done for different equalizer windows 
ranging from 1 to 7 taps. It can be stretched further. It has been found that when 
25 m = 8, four of the eight specified training sequence code (TSC) satisfy equation (11), and 
the auto-correlation matrix for the other four training sequences is in a specific sparse 
form: 
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16 0 
0 16 



U 



0 • • 16 • 
±4 0 • • 16 



(12) 



The inverse of equation (12) has the same sparse, symmetric structure: 



a 0 
0 1 



P 



a 



(13) 



where a = 1.0667, p = ± 0.2667. This means that if equation (9) is used for an estimation 
of 8-tap channels, only minor justification is necessary for certain TSCs: 



m ay 



(14) 



Property II indicates that not only a 5-tap channel, but all the 1-8 tap channels can 
10 be estimated efficiently using equation (9) without matrix inversion if any consecutive 16- 
symbol segment of the training sequence is used. Utilization of these properties has a 
fundamental impact on the computational complexity of the channel estimation for 
different equalizer window size. 

When an LSE approach is used for joint synchronization and equalizer window 
15 sizing, for a hypothetical sync position p and an equalizer window size m, by utilizing the 
properties discussed above, channel estimation can be computed as in equation (9). With 
the exception of the constant factor, it can be further expressed in an explicit convolution 
form, 

20 This can be considered as a FIR-style computation. For the next tap with k+ 1, a recurrent 
relation can be easily established, 
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15 

h(p,k + 1) « Z t(m+i-(k+l))r(p+m+0 
m ^ t(m+ i- k)r(p+ 1+ m+i) 

«= h(p + l.k) + t(m- 1 - k)r(p + m) - t(m + IS - k)r(p + m+ 16) 



Using this recursion in an HR-style computation, a new tap can be calculated using only 4 
real MAC operations (since the training sequences are real). This two-dimensional (sync- 
5 point tap-position) recursive relation can be depicted as shown in Fig. 5 for m = 5, w = 8. 
Computation of the last row and the last column is done by using equation (15), the rest 
can be recursively calculated using equation (16). 

To avoid channel estimation error caused by ISI corrupted receiving training 
symbols, by property I, we can always choose a training sequence segment that is close to 
10 the end. In practice, we often choose the last but one m + 15 symbols of the training 
sequence in the estimation (the final symbol in the training sequence is not used, to avoid 
possible precursor ISI corruption). When the next last w+15 symbols are chosen, equation 
(15) can be rewritten in a form independent of the equalizer window size m 

(17) h(p t k)=l t(i-k+9)r( P U + 9) 9 2<£<M 

i=0 max 

15 

Equation (16) can be adjusted accordingly by simply replacing m with 9. 

With the present invention, an efficient initial channel estimation method has been 
obtained for different equalizer window size spanning from 1-8 taps by utilizing the 
unique properties of GSMZEDGE training sequences described above. An efficient 
20 channel estimation at different sync points for different window size is essential to a 
practical implementation of joint optimization of synchronization and equalizer window 
sizing. The computational complexity of initial channel estimation by the approach 
according to the present invention is less than 10% of that of a brute force LSE approach. 
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The presented method has been put into simulation and verified. Simulation shows 
that a significant performance gain is observed in long time dispersive channels with joint 
optimization compared to the correlation based synchronization. 

It should be emphasized that the term "comprises/comprising" when used in this 
5 specification is taken to specify the presence of stated features, integers, steps or 
components but does not preclude the presence or addition of one or more other features, 
integers, steps, components or groups thereof. 

It should also be emphasized that while what has been described herein constitutes 
presently preferred embodiments of the invention, it should be recognized that the 
10 invention could take numerous other forms. Accordingly, it should be understood that the 
invention should be limited only insofar as is required by the scope of the following 
claims. 
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CLAIMS 

1. A method for estimating a transmission channel in a digital 
communications system which operates in accordance with GSM specifications, said 
communications system including a receiver which receives a signal transmitted over the 

5 transmission channel, said receiver including a channel estimation-based equalizer, said 
equalizer having a window size adapted to the actual delay spread of the transmission 
channel, and said equalizer using an estimate of the transmission channel for 
synchronizing and correcting said received signal, said received signal including a 26- 
symbol training sequence, said method comprising: 
10 using any consecutive 16-symbol segment of said 26-symbol training sequence to 

estimate, without matrix operation, the transmission channel for different equalizer 
window sizes at different synchronization points. 

2. The method according to Claim 1, where tap values for different 
15 synchronization points and tap positions are obtained in a two-dimensional recursive 

computation step for joint optimization of synchronization and window sizing of the said 
equalizer using Least Square Error (LSE) criterion. 

3. The method according to Claims 1 or 2, wherein said digital 
20 communications system operates in accordance with GSM/EDGE specifications. 

4. The method according to Claims 1 or 2, wherein said different equalizer 
window sizes comprises equalizer window spans of from 1-8 taps. 

25 5. The method according to Claims 1 or 2, wherein said using step comprises 

using delayed (shifted) segments in the training sequence for said estimation. 

6. The method according to Claim 2, wherein said tap values for different 
synch points are calculated in the 2-dimensional "north-west" direction. 

30 
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7. The method according to Claims 1 or 2, wherein the consecutive 16-symbol 
segment used is close to the end of the 26-symbol training sequence. 

8. The method according to Claim 5, wherein the consecutive 16-symbol 
5 segment used is the last but one m + 15 symbols, where m is the equalizer window size. 

9. The method according to Claim 1, wherein said using step comprises using 
an Least Square Error (LSE) approach for said estimation. 

10 10. The method according to Claim 2, wherein said using step comprises using 

an LSE approach for joint synchronization and equalizer window sizing. 

11. A method for estimating a transmission channel in a digital 
communications system which operates in accordance with GSMMEDGE specifications, 
15 said communications system including a receiver which receives a signal transmitted over 
the transmission channel, said receiver including a channel estimation-based equalizer 
which uses an estimate of the transmission channel for correcting said received signal, said 
signal including a 26-symbol training sequence, said method comprising: 



20 estimate, without matrix operation, the transmission channel for different equalizer 
window sizes of from 1 to 8 symbols long. 

12. The method according to Claim 1 1, wherein said using step comprises 
using delayed (shifted) segments of the training sequence for said estimation. 



13. The method according to Claim 12, wherein said using step comprises 
using a Least Square Error (LSE) approach for channel estimation. 



using any consecutivel6-symbol segment of said 26-symbol training sequence to 
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14. An apparatus for estimating a transmission channel in a digital 
communications system which operates in accordance with GSM/EDGE specifications, 
said system including a receiver which receives a signal transmitted over the transmission 

5 channel, said receiver including a channel estimation-based equalizer which uses an 
estimate of the transmission channel for correcting said received signal, said received 
signal including a 26-symbol training sequence, said equalizer using any consecutive 16- 
symbol segment of said 26-symbol training sequence for estimating, without matrix 
operation, the transmission channel for different equalizer window sizes of from 1 to 8 
10 symbols long. 

15. An apparatus for equalization of transmission channel in a digital 
communications system which operates in accordance with GSM/EDGE specifications, 
said apparatus including a receiver which receives a signal transmitted over the 

15 transmission channel, said receiver including a channel estimation-based equalizer which 
uses joint synchronization and equalizer window sizing to optimize the equalizer 
operation. 

20 
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